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PREFACE 


Progress  in  theoretical  research  has  been  continuing  toward  a  better 
understanding  of  collision  processes  of  various  kinds  involving  atoms, 
molecules  and  electrons.  High-speed  computers  and  accumulation  reservior 
of  numerical  techniques  as  well  as  computed  and  measured  data  have  contributed 
to  creating  a  progressively  favorable  condition,  from  which  new  ideas  are 
tested,  and  further  progress  is  made. 

In  Part  I  of  this  report,  we  have  applied  the  close-coupling  method  to 
the  electron-^  collision  problem  and  calculated  excitation  cross  sections 
of  the  Lymann-Birge-Hopfield  band.  Although  the  method  of  close-coupling 
remains  to  be  a  very  powerful  theoretical  tool,  its  effectiveness  can  be 
further  increased  substantially,  if  its  advantageous  features  are  combined 
with  those  of  the  R-Matrix  method.  This  is  described  in  Part  II  with  the 
electron-H2  collision  as  an  example. 

In  Part  III  attention  is  directed  to  slow-electron  collision  with  molecule, 
in  which  the  characteristic  feature  of  target  polarization  is  treated  by 
means  of  the  polarized-orbital  method.  Part  IV  deals  with  the  transition 
probability  of  the  dipole-forbidden  atonic  oxygen  transition  1Sq  -+  *D2 
(5577  A) ,  which  is  an  important  parameter  in  many  applications  in  the 
earth's  upper  atmosphere.  The  transition  probability  has  been  calculated 
with  multiconfiguration  Hartree-Fock  wave  functions. 
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PART  I 

ELECTRON  -  IMPACT  EXCITATION  CROSS  SECTIONS  OF  THE 

OF  THE  N7  MOLECULE  BY  THE  CLOSE -COUPLING  METHOD 

g  2 


1 
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I .  INTRODUCTION 

For  a  variety  of  research  fields  (e.g.,  geophysics,  laser  physics), 
the  electron -impact  excitation  of  electronic  states  of  molecules  provides 
basic  information.1"3  Theoretical  methods  of  calculating  excitation  cross 

sections  for  electronic  states  of  diatomic  molecules  range  from  the  Bom 

4  5  6  7 

approximation,  ,b  the  impact -parameter,  and  distorted- wave  methods  to 

8  9 

the  close-coupling  method.  ’  This  paper  reports  excitation  cross 

sections  of  ^(X1?:*  -*•  a1!!  )  by  a  two- state  close-coupling  (CC)  calculation  for 
2  2 

12.5  -  50  eV  range  of  incident  electron.  Within  the  framework  of  a  two- 
state  close-coupling,  the  calculations  are  performed  without  further 

9-12 

approximations.  In  particular,  although  there  are  well-reasoned  techniques 
of  approximating  electron  exchange,  no  recourse  is  made  to  them.  The  cross 
sections  are  averaged  over  the  orientation  of  the  molecule  with  respect  to 
the  direction  of  incident  electron.  This  is  equivalent  to  averaging  the  cross 
sections  over  the  initial  rotational  substates  and  summing  over  the  final 
rotational  substates,  to  a  very  good  approximation.13  The  electronic  states 
are  assumed  to  be  "vibrationless"  with  the  vertical  excitation  energy 
AF  =  9.1  eV.  Thus,  the  cross  sections  to  a  particular  vibrational  level  of 
the  a1^  state  may  be  computed  by  using  the  Franck -Condon- factor  approximation. 

The  formulation  of  the  close -coupling  method1^  is  briefly  recounted  in 
Sec.  II,  and  the  details  of  calculation  is  presented  in  Sec.  III.  In  Sec.  IV 
comparisons  are  made  with  the  Bom- type  calculations  and  with  experimental 
cross  sections.  A  short  discussion  in  Sec.  V  concludes  this  paper. 


I I .  CLOSE-COUPLING  METHOD 

Since  a  detailed  account  of  the  close-coupling  method  as  applied  to 
electronic  excitation  of  molecules  is  already  presented  in  Ref.  8  (hereafter 
referred  to  as  CL),  only  a  brief  summary  is  outlined  in  this  section.  As 
in  CL  the  molecule -fixed  frame  of  reference  is  used  here.  Then  the  scattered 
electron  and  target,  taken  together  as  a  collision  system,  are  characterized 
by  the  component  of  the  angular  momentum  along  the  molecular  axis  A  and  the  spin 
quantum  numbers  S,M.  The  target  wave  functions  are  appropriate  linear  com¬ 
binations  of  the  Slater  determinants  consisting  of  one-electron  orbital  functions 
*^1?),  where  X^  is  the  projection  of  angular  momentum  on  the  molecular 
axis  and  n^  is  used  to  distinguish  orbitals  of  same  symmetry  A^ .  The 
scattered-electron  function  is  decomposed  into  partial  waves  characterized 
by  angular  momenta  i  and  m  (=  A  -  X^).  A  short  notation  y  =  (nXX.)  is  used  to 
designate  a  channel. 

In  the  notation  of  CL  the  integro -differential  equation  for  the 

scattered -wave  function  F  ,  is: 

UP 


r 


dr 


=  2Z  ,,  [U  ,  „(r)  +  W  ,  ,,(r)]  F  „  (r) 


1 


where  U  and  W  are  respectively  the  direct  and  exchange  potentials  as  defined 
in  Eqs.  (16)  and  (17)  of  CL.  The  second  term  on  the  right-hand  side  of  Eq.  (1) 
is  to  ensure  that  F  ,  are  orthogonal  to  the  target  orbitals  by  means  of 
the  Lagrange  multipliers  A  .  Here,  the  two-center  molecular  orbital 

j  j 

functions  are  expanded  about  the  center  of  the  molecules  as 
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q  (nU  -  nU'l')  «  EA  QA(nAX,  ->  n’X’l'),  (6) 

and 

QCnA  v  n'X')  =  +x  f  t  QA(nX)l  ■>  n'A'A') .  (7) 

A—  8-=|  A-X  t  £,'  =  |  /\- A'  | 

In  order  to  facilitate  the  computation  in  practical  applications,  the 
spherical  coordinate  system  is  adopted  in  the  treatment  of  electronic 
excitation  of  molecules, inspite  of  the  fact  that  diatomic  molecules  have  only 
an  axial  symmetry  rather  than  a  spherical  one.  Thus,  in  this  approach,  a 
truncation  of  partial  waves  l  is  inevitable.  The  simulations  over  &  and  V 
in  f/i .  (7)  underscore  this  point.  This  trunication  can  be  recovered  with  the 
aid  of  a  partial -wave  analysis  of  the  Born  approximation.  For  a  given 


electronic  excitation  (nX  -*•  n'X'),  the  partial  cross  sections  by  CC,Q*-C(^ h{l,V 
and  by  Bom  approximation, A(£,£‘)  differ  substantially  for  small  (££'), 
but  approach  each  other  for  large  {IV) .  This  permits  one  to  compute  partial 
cross  sections  for  (££')  up  to  a  suitably  large  L,  and  use  the  Bom  partial 
cross  sections  for  (££')  >  L,  with  the  result 


»(CC) 


L 

=  EL 

u'  ' 


,  (CC) A 


+£^'>L 


£aQ(B)V,£') 


-  Q(BT)  +  -  QCB)iW,)l) ,  (8) 

where  1  is  the  total  cross  section  by  the  Bom  approximation.  This 

procedure  is  found  to  be  quite  satisfactory  here  as  it  was  in  other 

8,15 


e-molecule  collision  problems. 
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III.  DESCRIPTION  OF  CALCULATION 


A.  N7  wave  functions 

The  ground  (X^X+)  and  excited  fa^li  )  states  are  derived  from  the 
!,  g 

fol lowing  configurations, ^ J 

xV:  (la  )2(la  )2(2a  )2(2a  )2(1tt  )4(3a  )2, 
g  g  U  g  u  u  g  ’ 


a1!!  :  (la  )2(la  )2(2a  )2(2o  )2(1tt  )4(3o  )  (Itt  )  . 
g  g  uJ  g1  u'  u'  g'1  gJ 

17 

The  self-consistent- field  (SCF)  method  is  used  to  calculate  the  wave 

functions  by  applying  it  to  the  X2£+  state  and  to  the  a2n  state  at  the 

equilibrium  separation  (Rq  =  2.0675  aQ)  of  the  ground  state.  For  the 

basis  functions,  we  have  used  the  set  of  Gaussian-type  orbitals  (GTO) 

consisting  of  ten  ls-type  and  six  2p-type  GTO's  which  are  contracted  to 

18 

five  and  four  contracted  GTO's  respectively  as  given  by  Dunning.  In  the 

notation  of  Ref.  18,  this  set  reads  (10,6)/[5,4] .  The  total  energy  of 

-108.90376  a.u.  is  obtained  for  the  X2E+  state,  and  -108.43429  a.u.  for 

g 

the  a^rtg  state  in  the  present  work.  In  a  comparable  calculation  by  using  a 

19  20 

similar  (9,5)/[5,3]  GTO  basis  set  of  Dunning,  Truhlar  et  al.  obtained 

the  energy  (X^E+)  of  -108.8890  a.u.  They2^  improved  the  wave  function  by 
g 

including  two  3d-type  GTO's  in  the  basis,  i.e.,  by  using  (9  ,5  ,2) / f 4  ,5,2] 

set,  with  the  resulting  energy  of  -108.9732,  which  is  compatible  with 

21 

Nesbet's  earlier  calculation  with  the  Slater-type  orbitals  (STO)  as 
basis  (including  3d- type) .  The  qualities  of  these  and  other  wave  functions 
arc  compared  in  some  detail  by  Fades  et  al .  The  present  (10,6)/[5,4]  wave 
functions  are  expected  to  be  of  similar  quality  as  the  (9,  5) / [ 5 , 3]  functions 

4 

of  Ref.  20.  In  our  previous  Bom -approximation  calculation,  the  use 


7 


21  1 

of  Nesbet's  wave  functions  yielded  the  a  n  cross  sections  which  are 
about  6%  larger  than  the  results  obtained  by  using  the  CTO  functions 
similar  to  the  ones  used  in  this  work. 

Strictly  speaking,  the  orbitals  in  the  X'*'E+  function  do  not  satisfy  the 

s 

orthonormality  relation  with  those  in  the  a^II  function.  However,  each 

orbital  in  one  set  is  very  close  to  the  corresponding  one  in  the  other; 

in  the  worst  case,  the  overlap  of  3o  orbitals  of  the  X'*‘Z+  and  a^n  states 

is  C.9986.  The  orthogonality  between  orbitals  belonging  to  the  two  states 

-4 

is  found  to  be  within  3  x  10  except  for  the  overlap  0.0440  between  the 

2a  and  3a  orbitals.  Accepting  these  small  deviations,  we  used  all  the 
g  g 

orbitals  of  the  X^E+  state  as  common  to  both  states  except  Itt  orbital  in 
g  g 

order  to  facilitate  the  computation.  The  Gaussian  exponents  and  the 
corresponding  coefficients  of  the  orbitals  used  in  this  work  are  shown 
in  Table  I. 

B.  Coupling  Potentials 

The  method  of  computing  the  scattering  potentials  are  minutely  detailed 
in  CL;  in  this  subsection  the  particulars  pertaining  to  NL  (X^£+  ■+  a  n  )  are 

^  £  s 

discussed.  The  direct  coupling  potentials  are  due  to  interaction  of  scattered 
electron  with  the  nuclei  (V^) ,  and  with  the  bound  electrons  (Ve) ,  viz., 


V"-  -'(r)  '  '  V),  YM-»Cr,V.A-X>* 


UfM 


X  [|fA  -  ff1  *  |fB  -  ff1! 


(9) 


=  /  dr  V-A(r)  Yr,A-A'* 


Xjfj  /  <f)j  0"')  |r  -  r'  (r')dr' , 


x 


(10) 
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where  Z  (=7)  is  the  nuclear  charge,  and  r^,  Tg  are  the  position  vectors  of  the 
two  nuclei.  In  Bq.  (10)  the  summation  over  j  indicates  the  appropriate  pairs 

I 

of  molecular  orbitals  4k  and  «Jk  ,  and  f ^  are  purely  numerical  factors. 

These  potentials  are  computed  numerically  by  the  method  described  in  Secs. 

Ill  and  IV  of  CL.  All  14  electrons  of  ^  are  fully  taken  into  account  in  this 
work. 


lEq. 


The  exchange  potentials  are  integral  operators,  which  may  be  written  as 
(46)  of  CL] 

V(r,Fpv(r)  = '  hgj  K?ghjm)r  /drY >c,gw 

x  Y*,A-\(r)h  (f)  l  rr'\.u.,(r")<ir' 

oo  _ 


+  rK  /  r'  K  1F  ,  „(r')dr')  r’ 
r  mm 


/  dr-  YKjg(r-)Y£t>A.A,(r-) 


x  0*  (?’)], 


(ID 


where  the  summation  over  j  encompasses  the  appropriate  pairs  of  orbitals 


3 


and  <Pj  with  numerical  factors  g^ .  Here,  the  range  of  K  is  unbounded.  This 
can  be  illustrated  by  considering 


l  /dr  YK,g  (r)  YM-ACrJ^j  (r) 

=  LJ 
K 


™  J  dr  Y  (r)  Y  (r)  Y  (r)4>'  (r) 

i  K,g  t,s’A  Vi  J-h 


(12) 


t 

where  the  molecular  orbital  4k  is  decomposed  as  shown  in  P.q.  (2).  Therefore, 
Bq.  (11)  can  be  evaluated  only  approximately,  i.e.,  with  a  truncation  in  the 
K- series.  The  results  presented  in  this  paper  are  based  on  retaining  two 
leading  terms  of  K  in  Bq.  (11).  However,  we  carefully  monitored  the  accuracy 
of  this  scheme  by  repeating  the  calculation  with  five  such  terms  included 


at  a  few  points  with  respect  to  both  the  incident -electron  energy  and  A. 

The  difference  in  cross  sections  between  the  two  sets  of  calculations  is 
quite  negligible  (M3. 3%),  indicating  that  the  two  leading  terms  of  K  in 
Eq.  (11)  were  sufficient.  As  with  the  direct  potential,  all  14  electrons  are 
fully  taken  into  account  in  the  exchange  potential. 

C.  Scattering  equation 

The  solution  of  Eq.  (1)  by  straightforward  numerical  integration  is 

not  possible  because  of  the  exchange  term  which  requires  a  prior  knowledge 

of  the  solution  F(r)  for  all  r  as  shown  in  Eq.  (11).  In  the  present  iterative 

procedure,  (1)  the  exchange  term  WF  of  Eq.  (11)  is  constructed  using  the 

F  generated  in  the  previous  cycle  of  iteration  (we  set  WF  =  0  for  the  first 

cycle);  (2)  using  the  WF  so  obtained,  new  solutions  F  are  calculated; 

(3)  steps  (1)  and  (2)  are  repeated  until  a  self-consistency  is  achieved 

between  the  successive  cycles.  This  procedure  works  well  for  large  A 

value  and/or  at  high  incident -electron  energy,  and  the  cross  sections  are 

converged  to  0.015  within  4-8  iterations  depending  on  A  and  the  incident 

energy.  However,  for  small  A  or  low  incident  energy,  this  method  fails 

23 

to  give  convergent  results.  In  such  cases  we  obtained  the  solution  by 
the  noniterative  integral  equation  method  (NIEM)  of  Smith  and  Henry.24  However, 
NIEM  is  a  much  more  costly  process  than  an  iterative  one  in  terms  of  computer 
usage.  Therefore,  we  used  the  iterative  scheme  whenever  possible  (usually  for 
A  >  3) ,  and  for  other  cases  (A  =  0  -  3)  we  used  NIEM.  The  step-sizes  (in  units 
of  aQ)  used  in  the  numerical  integration  are  6r  =  0.0125  for  0  <  r  <  1.5, 

6r  =  0.025  for  1.5  <  r  <  3.0,  <5r  =  0.05  for  3.0  <  r  <  5.0,  and  6r  =  0.1  for 
5.0  <  r  <  100.  The  effective  range  of  exchange  potential  is  about  the  radial 
extent  of  molecular  orbitals.  In  this  work  the  exchange  potentials  are 
retained  to  r  =  8.0  aQ,  beyond  which  they  arc  totally  negligible.  The  scattering 


10 


matrix  S  ,  in  Eq.  (3)  is  determined  at  three  different  points  (r  =  SO,  75, 
100  aQ) ,  and  the  cross  sections  are  obtained  at  each  point.  No  significant 
change  in  cross  sections  is  found  with  respect  to  the  "boundary-matching" 
points . 


1 

1 
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IV.  RESULTS 

The  excited  state  a^II  is  a  degenerate  state  with  \  =  ±1.  We  made  two 

o 

separate  two -state  calculations  including  one  degenerate  state  each  with  the 
ground  state.  The  other  option  wuld  be  to  make  a  three-state  calculation. 

We  expect  similar  results  from  either  set  of  calculation.  In  an  analogous 

X  +  3  8 

case  of  FL  (XT  ■*  c II  )  the  difference  was  M)%. 

2  g  uJ 

Since  is  a  homonuclear  molecule,  the  scattering  equations  separate 
into  two  blocks  of  even  (gerade)  and  odd  (ungerade)  symmetry.  The  partial 
cross  sections,  as  defined  in  Eq.  (5) ,  at  the  incident-electron  energy 
E  =  50  eV  are  shown  in  Table  II.  For  a  singlet-singlet  excitation  the  partial 
cross  sections  (in  series  of  A)  converges  slowly  as  seen  from  Table  II.  To 
achieve  convergence,  the  scope  of  calculation  could  be  expanded,  in  principle, 
with  respect  to  £  and  A.  An  alternative  to  this  procedure  is  to  augment 
the  close -coupling  calculation  with  the  Born  approximation  (BA)  as  described 
in  Sec.  II.  The  rationale  of  this  procedure  is  based  on  the  observation  that 
the  high  £  partial  waves  cannot  penetrate  close  to  the  target  so  that  the 
effect  of  electron-exchange  is  diminished  to  the  point  of  insignificance. 

Under  this  condition  the  BA  is  expected  to  become  valid.  This  is  illustrated 
in  Table  III,  where  comparisons  are  made  of  partial  cross  sections  computed  by 
CC  and  by  BA  at  E  =  50  eV.  It  is  seen  that  the  agreement  is  quite  good  beyond 
£  >  7.  In  other  words,  the  objective  of  CC  is  to  calculate  low  £  partial 
cross  sections,  which  the  BA  overestimates,  as  shown  in  Eq.  (8).  The 
quantities  appearing  in  Eq.  (8)  are  shown  in  Table  IV  at  five  different  energies. 
The  second  column  shows  the  total  Born  cross  sections.  From  the  Born  total, 
we  replace  the  Born  partial  (third  column)  by  the  corresponding  CC  partial 
(fourth  column)  to  obtain  the  cross  sections  (last  column).  As  mentioned 
before,  we  included  in  the  calculation  interactions  between  the  incident 


and  all  14  target  electrons  of  ^  molecule.  However,  with  a  view  toward 

future  works,  we  performed  test  calculations  at  E  =  17  eV,  A  =  1,  and  even 

symmetry  with  (a)  la  and  la  omitted,  and  (b)  lo  ,  lo  ,  2o  ,  and  2c  omitted. 

'  1  g  u  g  u  g  u 

2  2 

The  resulting  cross  sections  are  (a)  0.0472  aQ  ,  and  (b)  0.0587  as  compared 

2 

with  the  full  calculation  value  of  0.0487  aQ  .  It  is  seen  that  the  tightly- 
bound  la  and  lou  orbitals  could  have  been  excluded  from  the  calculation  without 
significantly  03%)  changing  the  final  outcome. 

The  cross  sections  in  Table  IV  are  also  plotted  in  Fig.  1,  and  compared 

2  cj  26 

with  the  Born-Ochkur  results,  and  experimental  measurements.  ’  The 
Born-Ochkur  approximation  (B0)  yields  good  results  down  to  ^30  eV.  Since  the 
success  of  BO  is  largely  phenomenological,  it  could  be  used  only  for  a  rough 
estimates  of  cross  sections. 

The  X1?/  -►  a1!!  excitation  has  been  studied  experimentally  by  numerous 

o  o 

25-3? 

groups.  ‘  In  Fig.  1  we  include  two  of  the  latest  results  by  Cartwright 

25  2^)  27  32 

et  al.  and  by  Finn  and  Doering.  The  earlier  experiments  L  are  discussed 

25  26 

in  some  detail  in  the  these  papers.  ’  There  are  some  differences  both  in 

magnitude  and  shape  of  the  excitation  function  between  the  two  sets  of  experimental 

25 

data  shown  in  Fig.  1.  Cartwright  et_  al/  attribute  the  difference  in  magnitude 
to  the  cascade  contributions  which  were  not  considered  by  Finn  and  DoeringV 
The  other  experimental  excitation  functions  (also  without  cascade  correction) 
lie  generally  between  the  two  curves  shown  in  Fig.  1. 

The  present  theory  agrees  reasonably  well  with  experiments  above  30  eV, 
but  the  disagreement  below  30  eV  is  apparent.  Although  there  are  considerable 
uncertainties  in  experiments,  all  experimental  excitation  functions  show  the 
peak  occurring  somewhere  between  15  and  18  eV.  Therefore, 
a  good  part  of  the  discrepancy  between  theory  and  experiment  below  30  eV 
is  likely  due  to  the  shortcomings  of  the  two-state  close-coupling  theory. 
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From  the  standpoint  of  the  theoretical  formalism  employed  in  a  calculation, 

the  two-state  close -coup ling  is  the  most  comprehensive  one  applied  to 

calculate  the  cross  sections  of  the  a1!!  state  to  date,  and  the  Bom, 

g 

Bom-Ochkur,  or  distorted-wave  method  cannot  be  expected  to  outperform 

the  close -coup ling,  since  these  methods  involve  further  approximation  to 

the  two-state  c lose- coupling.  Improvement  of  wave  functions  is  always 

desirable;  it  is  more  difficult  to  answer  to  what  extent  the  computed 

cross  sections  are  affected  by  the  improvement  of  wave  functions.  For 

an  estimate,  we  recall  that  the  inclusion  of  3d-type  function  affected  the 

4 

cross  sections  by  6%  in  the  Bom-Ochkur  calculation.  The  configuration 
interaction  (Cl),  as  a  means  of  accounting  for  electron  correlation,  also 
needs  to  be  considered.  However,  Billingsley  and  Krauss33  report  that 


I  +  _  2^ 

the  value  of  quadripole  moment  of  N2(X  £g)  changes  from  -1.29  x  10 

_  2 

(single  configuration  function)  to  -1.22  x  10  esu-cm  (18  configuration 

function) ,  which  indicates  much  less  influence  of  Cl  in  N2  than  in  02 

where  the  oscillator  strength  of  the  Schumann- Runge  system  is  affected 
34 

by  a  factor  of  3.  Therefore,  we  believe  that  the  improvement  in  the 
wave  functions  (Cl  included)  will  not  change  the  computed  cross  sections 
by  the  two- state  close-coupling  to  the  extent  of  bringing  theory  and 
experiment  to  a  satisfactory  agreement. 
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V.  CONCLUSION 

In  sunmary,  two-state  close -coupling  calculations  have  been  performed 

to  obtain  the  excitation  cross  sections  of  a^n  state,  augmented  by  the  Bom 

S 

partial  wave  analysis  to  handle  high  Si  partial  waves.  The  present  cal¬ 
culations  and  experiments  are  in  good  accord  above  30  eV.  In  the  experi¬ 
mental  peak  region  the  weight  of  experimental  evidence  points  to  under¬ 
estimating  of  cross  sections  by  the  two-state  close -coupling.  In  view 
of  the  discussion  presented  in  the  last  paragraph  of  Sec.  IV,  the  target 
polarization  (i.e.,  the  distortion  of  ^  as  the  incident  electron  approaches 
it)  emerges  as  an  important  factor  to  be  considered  toward  a  better 

description  of  the  X^E+  -*•  a^n  excitation  process. 

S  § 

To  incorporate  the  target  polarization,  by  such  methods  as  that  of  the 
35  36 

polarized  orbitals,  *  °  would  add  another  dimension  of  complexity  to 
the  present  theoretical  treatment  of  electronic  excitation  of  molecules. 
However,  efforts  should  be  made  in  this  direction,  in  view  of  the  significant 


influence  of  the  polarization  effects  on  the  electron-atom  collisions. 


35-37 
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Table  1.  Gaussian  exponents  and  coefficients  of  Molecular  orbitals 
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Table  II.  Partial  cross 


in  units  of  0.1 


A 

aV 

g 

even 

(  A  =  +1) 

odd 

0 

0.230831 

0.205208 

1 

0. 146895 

0.137790 

2 

0.214090 

0.100259 

3 

0.134843 

0.207456 

4 

0.130350 

0.119458 

5 

0.066402 

0.077523 

6 

0.033294 

0.040752 

7 

0.013894 

0.019317 

8 


0.002889 
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tions  as  defined  in  Eq.  (5) 

2  at  E  =  50  eV. 
o 


aV 

g 

even 

(  A  =  -1) 

odd 

0.105154 

0.154361 

0.094702 

0.110489 

0.067513 

0.048896 

0.023791 

0.035915 

0.016711 

0.012731 

0.004881 

0.008817 

0.002909 


0.002669 
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Table  III.  Partial  cross  sections  as  defined  in  Eq.  (6) 

2 

in  units  of  0.1  a  at  E  =  50  eV. 
o 


Partial  wave 
(incident) 

Q 

cc 

(ft,  ft*  =  t-2) 

Born 

Q  (M 

cc 

=  A) 

Bom 

1 

0.107902 

0.262280 
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0.007377 

0.013718 

0.105873 
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0.043471 

0.063359 
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0.290970 
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0.220897 
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0.131745 

0.139521 

7 

0.210181 

0.220489 

0.082937 

0.085231 

8 

0.167138 

0.169191 

0.052182 

0.052686 

9 

0.127740 

0.128287 

0.033137 

0.033227 
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Table  IV.  Total  and  partial  cross  sections  in  units 
2 

of  aQ  by  Bom  and  CC  as  indicated  in  Eq.  (8). 


Energy  (eV) 

Bom  (Total) 
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CC  (Partial) 
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IIG.  1.  Excitation  cross  sections  of  the  a  ng  state  calculated  by 
CC  (solid  line);  by  B0  (long-short  dashed  line) ;  and  experimental 
cross  sections  of  Ref.  25  (uniform  dashed  line);  and  of  Ref.  26 
(long-short-short  dashed  line). 
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I .  INTRODUCTION 

For  low-energy  electron  scattering  the  exchange  interaction  between  the 
incident  electron  and  the  target  electrons  plays  an  important  role.  Particu¬ 
larly  in  the  case  of  excitation  of  the  target  atom  (or  molecule)  with  a 
change  in  the  electron-spin  multiplicity,  electron  exchange  is  the  major 
interaction  responsible  for  this  type  of  process.  While  the  exchange  inter¬ 
action  has  been  treated  in  an  ab  initio  manner  in  the  theory  of  electron-atom  and 
electron-molecule  excitation  based  on  the  method  of  close  coupling  (CC),*  5  the 

computing  cost  for  this  type  of  calculation  is  usually  quite  large.  There 
4 

have  been  attempts  to  approximate  the  electron -exchange;  however,  the  degree 
of  success  of  these  models  must  be  judged  ultimately  by  comparison  with  the 
results  of  the  exact  treatment  of  electron-exchange. 

The  introduction  of  the  R-Matrix  (RM)  method5  to  atomic  physics  is 
largely  based  on  the  observation  that  the  effective  range  of  the  electron- 
exchange  is  confined  to  the  radial  extent  of  the  target  orbitals  (say,  r  =  a) . 

The  R-Matrix  radius  a  is  chosen  so  that  beyond  r  =  a  the  electron -exchange  may 
be  completely  ignored.  Inside  the  boundary  the  electron-exchange  is  treated 
rigorously.  However,  by  limiting  the  radial  extent  (r  <_  a) ,  it  becomes 
possible  to  represent  the  continuum  functions  by  a  finite  number  of  basis 
functions.  Then  the  scattering  equations  within  the  boundary  can  be  solved 
more  efficiently  by  using  the  matrix  method  in  contrast  to  having  to  solve 
coupled  integro-differential  equations  in  CC.  Outside  the  sphere  the  exchange 
terms  are  neglected  so  that  the  CC  scattering  equations  involve  no  integral 
operators  and  can  be  handled  easily. 

In  the  practical  application  of  RM,  the  choice  of  basis  functions 
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becomes  an  important  factor  with  regard  to  both  the  efficiency  of  computation 

and  the  accuracy  of  the  results.  Various  techniques  have  been  used  to  generate 

basis  functions,  including  analytic  floating  Gaussian  bases^  and  approximate 

5  7 

solutions  to  the  static  potential  problem  without  exchange.  *  This  report 
describes  a  practical  procedure  by  which  approximate  CC  solutions  are  used 
as  a  basis  for  the  RM  calculation.  Such  a  hybrid  CC-RM  procedure  affords  the 
advantage  that  a  very  small  basis  provides  an  excellent  representation  of 
the  scattering  solutions  within  the  R-Matrix  boundary,  and  hence  greatly 
reduces  the  cost  of  the  calculation. 

Secs.  II  and  III  will  outline  same  key  features  of  CC  and  RM  respectively. 
No  attempt  will  be  made  to  provide  a  complete  derivation  of  either  procedure, 
as  such  derivations  can  be  found  in  the  references  cited.  Our  emphasis  will 
rather  be  on  Sec.  IV,  in  which  the  scheme  for  combining  the  two  procedures 
will  be  described.  In  Sec.  V  a  comparison  is  made  of  the  two  sets  of 
excitation  cross  sections  of  the  b3£+  state  of  H-  computed  by  the  hybrid 

U  t 

CC-RM  procedure  and  by  the  standard  CC  method. 


.  J  •„ ...j.  .  .. 
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II.  CLOSE-COUPLING  METHOD 


The  details  of  the  close-coupling  method  as  applied  to  electron-molecule 
scattering  have  been  presented  in  Ref.  2,  and  will  not  be  repeated  here.  In 
the  molecule -fixed  frame  of  reference,  the  scattered  electron  and  target, 
taken  together  as  a  collision  system,  are  characterized  by  the  component 
of  the  angular  momentum  along  the  molecular  axis  A  and  the  spin  quantum 
numbers  S,M.  The  target  wave  functions  are  appropriate  linear  combinations 
of  the  Slater  determinants  consisting  of  one-electron  orbital  functions 
4>j(n-A.|r)  where  A  j  is  the  projection  of  angular  momentum  on  the  molecular 


axis  and  n^  is  used  to  distinguish  orbitals  of  the  same  symmetry  A . .  The 
scattered -electron  function  is  decomposed  into  partial  waves  of  angular 
momenta  sl  and  m  (=  A  -  Aj).  The  shorthand  notation  y  =  (nA£)  will  be  used  to 
designate  a  channel.  The  set  of  integro-differential  equations  for  the 
scattered  wave  functions,  F  .  are 

y'y 


d?  '  *  V'1J  *  k'2)  F‘‘'»(r)  '  2  V[UpV(r)  *  WuV'(r)1P1VuM 


+  (orthogonality),  (1) 

where  U  and  W  are  the  direct  and  exchange  potentials  respectively,  and 
"(orthogonality)"  represents  the  constraints  included  to  ensure  that  the  F 

p’l 

are  orthogonal  to  the  target  orbitals.  The  direct  coupling  potentials  are  due 
to  the  interaction  of  the  scattered  electron  with  the  nuclei  (l^) ,  and  with  the 
bound  electrons  (Ue),  i.e.. 


U  ,  (r)  =  ,  +  ue 

MM’  mm'  yy'  * 


(2) 


u„y'  W  =  "6(nA) ,  (n'A ')  Z  /  dr  Y^.^r)  V,A-A’(?) 

x  [I?.-?!'1  + 


(3) 
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B 


■  /*YM-»(r)  Va-X'W 


*  Zjfj  ^  ')  I  r ' -r | -1  ♦jtr'Jdr',  (4) 

where  Z  is  the  nuclear  charge,  and  r^,  are  the  position  vectors  of  the  two 
nuclei.  In  Eq.  (4)  the  summation  over  j  indicates  the  appropriate  pairs 
of  molecular  orbitals  (MO)  4> ^  and  <t> j  ,  and  f \  are  the  corresponding  numerical 
factors  (see  Ref.  2  for  more  details) .  The  exchange  potentials  are  integral 
operators ,  viz . , 


V(r)FuV*  (r)  =  '  Zjgj  ZK,g  iTOfl  r/d?  YK,g^Y*,A-X^^j(?) 

*  ([r'K'1/Jr'KFiJ,lill(r')dr'  -  rK/Jr,'K'1Fyly„(r’)dr' 

+  rKj«r,-K-lF^^(r,)drt]  r,jd?,  YK>g(r»)Y VtA.V  (5) 


where ,  as  in  Eq.  (4),  the  summation  over  j  covers  the  appropriate  MO's  <Jk 

and  with  the  nunerical  factors  g^  .  Because  of  the  third  integral  in  the 

square  bracket,  the  solution  of  the  integro-differential  equations  in  Eq.  (1) 

2  3 

takes  a  great  deal  of  computing  time.  ’  This  is  a  discouraging  aspect  of 
CC  in  practice. 


•3* 


X. 
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III.  R-MATRIX  METHOD 
5 

Burke  and  coworkers  introduced  the  R-Matrix  method  to  the  electron- 

atom  and  electron -molecule  scattering  problems  with  detailed  derivation  of 

the  formulation.  It  is  motivated  by  the  observation  that  the  exchange  terms 

have  limited  ranges  as  evidenced  by  the  first  integral  in  Eq.  (5) .  The 

R-Matrix  radius  a  is  chosen  such  that  the  exchange  is  negligible  beyond 

r  =  a.  Thus,  outside  of  this  boundary  the  differential  equation  of  Eq.  (1) , 

with  WF  =  0,  can  be  solved  quite  expeditiously. 

Within  the  boundary,  the  scattered  waves  F  are  expanded  in  terms 

of  chosen  basis  functions  v  . 

pi 


-  .  (r)  =  l.  ,c  , . , 

p  ’  JJ  1  T  JJ  *  1  f 


,u  Vf(r)- 


(6) 


It  is  shown  in  Ref.  5  that  diagonalizing  the  Hamiltonian  matrix  leads  to  an 
equation 


-  ,  (a)  =  a  E  „R  ,  ,,F'  ,,  (a), 

p  p  p  PM  p  p 


I?) 


where  F'(a)  is  the  derivative  of  F(r)  at  a.  The  entire  interaction  of  both  the 
direct  and  exchange  types  is  absorbed  in  Eq.  (7),  which  relates  the  scattered- 
wave  functions  and  their  derivatives  at  r  =  a.  The  R  matrix  is  shown  to  be'’ 


R  =  (l/2a)  vT(a)  (IIs  -  ES)"1v(a),  (8) 

where  If’  is  the  symmetric  part  of  the  Hamiltonian  matrix  and  E  is  the  diagonal 
energy  matrix.  The  overlap  matrix  S  arises  due  to  the  use  of  nonorthogonal 
basis  functions  v. 

Eq.  (7)  is  exact  and  is  limited  only  by  the  accuracy  of  the  R  matrix.  The 
latter  in  turn  depends  solely  on  the  basis  functions  v  employed,  in  Eq.  (8). 


.  Wlh._.L 


Thus,  we  see  that  it  is  extremely  important  to  devise  a  means  of  securing  a 

small  but  representative  basis  set  in  RM.  However,  this  aspect  is  most 

5  -7 

wanting  in  the  applications  of  RM,  although  several  schemes  have  been 
employed.  In  the  following  section,  we  combine  the  advantageous  features 
of  the  RM  and  CC ,  thereby  gaining  a  substantial  improvement  in  efficiency 
over  the  conventional  RM. 
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IV.  HYBRID  CC-RM 

In  the  conventional  CC  approach,  the  solution  of  the  integro-different ial 
CC  equations  is  quite  time-consuming  in  the  vicinity  of  the  molecule  where 
the  exchange  interaction  must  be  included.  This  region  is  exactly  the  region 
RM  is  designed  to  handle.  The  most  difficult  aspect  of  RM  is  to  generate  a 
small  but  adequate  basis  set.  On  the  other  hand  CC  is  well  suited  for  providing 
such  a  basis  set  by  approximate  solution  of  the  CC  equations.  Thus  a 
combination  of  the  CC  and  RM  may  lead  to  significant  improvement  over  each 
individual  method.  Our  prescription  for  a  hybrid  CC-RM  procedure  is  as  follows: 

(1)  The  integro-differential  equations  [Eq.  (1)]  of  CC  are  solved 
iteratively  for  two  iterations  as  described  in  IIIC  of  Ref.  3. 

(2)  The  resulting  approximate  solutions  denotes  approximation 

to  F  ,  )  are  used  as  basis  functions  v  to  calculate  the  R  matrix  in 

PM 

Eq.  (8). 

(3)  The  differential  equations  for  r  >  a  are  solved  without  exchange  with 
Eq-  (7)  providing  the  starting  solutions. 

We  describe  below  each  of  the  steps  in  same  detail. 

A.  Generation  of  basis  set 

In  the  conventional  CC  the  iterative  solutions  of  Eq.  (1)  take  four  or 
more  iterations  to  converge  under  favorable  conditions.^  However,  in  the 
hybrid  CC-RM  approach  we  seek  only  approximate  solutions  of  Eq.  (1)  and  use 
them  as  basis  functions,  thus  two  iterations  are  found  to  be  adequate.  Further, 
these  solutions  are  obtained  with  the  orthogonality  constraints  with  respect 
to  the  bound  MO's  so  that  the  computation  of  the  R  matrix  in  Eq.  (8)  is 
simplified.  The  basis  set  for  a  particular  exit  channel  u'  is  formed  from 

with  y  =  1,2,. . . ,n  corresponding  to  all  incident  channels.  These  functions 
4^,^  are  normalized  within  the  R-matrix  boundary  (0,a)  and  redundant  ones  are 
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discarded  by  an  overlap  criterion.  However,  the  dominant  diagonal  ones 
ij.  are  always  retained. 

w 

B.  Solutions  in  the  exterior  region 
The  RM  procedure  has  served  its  purpose  when  R(a)  has  been  obtained. 
The  CC  equations  without  exchange  are  to  be  solved  in  the  exterior  region 
from  r  =  a.  That  is,  for  r  >  a  Eq.  (1)  is  reduced  to 


+  k<2)  f  (r)  =  2  E  „  U  ,  ,,F  „  (r) . 
Ti  p  y  y  y  y  y  y 


(9) 


The  orthogonality  constraints,  like  the  exchange  terms,  can  be  ignored  here. 
8 

The  Numerov  method  is  used  to  generate  the  solutions  until  the  asymptotic 
form  has  been  reached ,  i . e . , 


Fyty(r)  *  (lA^)[6y,ue'l(k,r-£,7T/2)-ei(k'r'riT/2)S^].  (10) 

The  cross  sections  are  extracted  from  the  scattering  matrix  S^y  in  the  usual 
2 

manner  » 


QSM^(nAs£  -*■  n’A’s’fc’)  =  (tt/K2)  1 5yty  -  sj^|2, 


(ID 


Q(nAs  -+  n'A’s')  =  E  (^(nAsit  -►  n'A's'X.').  (12) 

For  the  purpose  of  discussion  we  also  define 


cma  sma 

Q^v(nAs  +  n’A’s')  =  l  (finAsit  +  n’A’s’fc’). 

fcfc' 


(13) 


The  manner  in  which  R(a)  is  used  to  construct  the  starting  solutions  is  now 

discussed.  Near  the  boundary  r  -  a  the  true  solutions  F  ,  can  he  written  as 

HP 


Fp.pW  -  N[Ju,M(r)  *  YuV  WK^./r)] , 


(14) 


with 


Ju'u(r)  "  VA*-r 


k 


(15) 
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Y . W  =  6y'MkMJ'r  y£..,(ky*r)' 


y’y 


(16) 


where  and  are  the  spherical  Bessel  functions  of  the  first  and  second  kinds 
respectively.  Although  the  KCr^a)  matrix  differs  from  its  asymptotic  value 
K(r-«“),  we  expect  that  it  varies  slowly  with  r  in  this  region  so  that 


K(a+6r)  =  K(a),  (17) 

the  step-size  6r  being  of  the  order  of  0.1  a^.  The  normalization  constant  N  in 
Eq.  (14)  is  imnaterial  and  set  to  1.  By  substituting  Eq.  (14)  into  Eq.  (7), 
we  obtain 

K(a)  =  -  (Y  -  aRY')'1  (J  -  aRJ'),  (18) 

where  all  matrices  are  evaluated  at  r  =  a,  and  the  prime  indicates  the  derivative 
with  respect  to  r.  Here,  we  have  ignored  K’ (a)  in  view  of  Eq.  (17).  Thus, 
with  the  K  matrix  as  given  by  Eq.  (18),  we  construct  Fyty  at  two  starting 
points,  viz.. 


F  ,  (a) 
y'yv 


J  ,  (a) 
y  yv 


YyI(J,(a)Kyly(a), 


and 


(19) 


F  ,  (a+<5r)  =  J  ,  (a+5r)  +  Y  ,  ,(a+6r)K  , 

u'u  u'y  y  y  y  y 


(a), 


(20) 


With  these  two  starting  points,  the  three-point  recursion  formula  of  Numerov 
may  be  used  to  perpetuate  the  solutions  to  any  desired  r. 


irtfto.ni.i  A  V 
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V.  RESULTS 


The  hybrid  CC-RM  procedure  has  been  used  to  calculate  the  excitation 

(X^Z*  -*■  b^E*)  cross  sections  of  H_  at  incident  electron  energies  of  13,  15,  20, 

6  u  “ 

and  30  eV.  Since  these  cross  sections  have  been  already  calculated  by  the 

2 

close -coupling  method,  a  good  test  of  the  present  procedure  can  be  made  by 

comparison.  The  wave  functions  and  various  potentials  employed  in  this  work 

2 

are  identical  to  those  used  in  the  CC  calculation.  At  each  incident  energy 

the  approximate  CC  parent  basis  functions  are  generated  in  a 

numerical  tabular  form  for  channel  angular  momentum  up  to  Z  =  5.  Redundancy 

within  this  parent  basis  set  is  removed  using  a  normalized  overlap  tolerance 

of  0.9.  rhe  resulting  basis  set  is  very  small,  with  only  two  or  three  basis 

functions  for  each  exit  channel  generally.  Table  I  shows  a  typical  number 

of  basis  functions  used  for  the  calculation  of  A  =  0  and  E  =  15  eV.  The 

R  matrix  is  then  calculated  with  the  boundary  a  =  20  aQ. 

Table  II  shows  the  convergence  of  the  partial -wave  expansions  for  a 

typical  (A,E)  block  as  defined  in  Eq.  (11)  with  S  =  M  =  h.  The  partial  cross 

sections  due  to  the  largest  Z  retained  are  seen  to  be  four  to  five  orders 

of  magnitude  smaller  than  the  largest  partial  Cross  sections  in  the  block. 

Since  fewer  partial  waves  participate  at  lower  energies,  this  convergence  at 

E  =  30  eV  as  shown  in  Table  II  should  guarantee  convergence  at  all  lower  energies. 

The  convergence  with  respect  to  the  A-expansion  is  shown  in  Table  III  for 

the  incident  energy  E  =  15  eV  (near  peak).  Por  the  excitation  of  homonuclear 

molecules,  the  partial  cross  sections  as  well  as  the  differential  equations 

separate  according  to  even  and  odd  parity.  The  partial  cross  sections  of 

A  =  0,1,2  are  shown  to  be  quite  sufficient.  The  agreement  between  the  present 

2 

hybrid  CC-RM  and  the  standard  CC  results,  shown  in  Table  IV  and  Figure  1, 


is  within  5%  at  worst.  This  excellent  agreement  demonstrates  the  adequancy 
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of  the  small  CC  basis  set.  Test  calculations  with  the  spherical  Bessel 

functions  as  basis  show  that  it  takes  90  Bessel  functions  to  achieve  a  similar 

2 

degree  of  agreement  with  the  previous  close-coupling  calculation  as  the 
present  calculation.  The  difficulty  of  achieving  convergence  with  the  Bessel 
functions  stems  from  the  fact  that  the  Bessel  functions  are  in  reality 
components  of  a  plane  wave,  whereas  the  true  scattered  functions  are  considerably 
phase-shifted  and  otherwise  modified  from  the  plane  wave  in  the  interior 
reg  ion . 

In  general  the  R  matrix  boundary  a  is  chosen  by  compromising  the  two 
conflicting  requirements:  (1)  the  larger  the  a,  the  less  of  the  exchange  is 
ignored;  (2)  the  smaller  the  a,  the  better  the  expansion  Eq.  (6)  becomes  with 
a  given  number  of  basis  functions.  The  rather  large  value  of  a  =  20  a^  is 
chosen  to  provide  a  very  stringent  test  of  the  ability  of  the  CC  basis  to 
adequately  represent  the  scattering  functions  in  the  interior  region.  The 
previous  CC  calculation,  to  which  we  compare  present  results,  also  included 
the  exchange  up  to  r  =  20  a^.  However,  the  boundary  a  can  be  significantly 
reduced.  For  example,  even  for  the  diffuse  la  orbital,  integration  beyond 
r  =  14  a^  does  not  contribute  more  than  10  ^  to  the  normalization  value  of 
unity.  Had  we  chosen  this  (14  a^)  or  yet  smaller  value  for  a,  the  rep¬ 
resentation  by  our  present  basis  functions  would  have  been  even  better.  This 
would  also  have  reduced  the  computing  cost  still  more. 


V.  CONCLUSION 


We  find  that  the  hybrid  CC-RM  described  in  this  paper  is  a  far  more 
efficient  method  than  either  CC  or  RM  method  used  separately.  From  the  point 
of  view  of  CC  this  method  eliminates  the  time-consuming  process  of  solving 
the  integro -differential  equations  in  the  interior  region  by  "borrowing"  the 
matrix  method  of  RM.  In  the  present  calculation  of  with  six  channels  or 
less,  the  overall  computer  usage  was  reduced  by  a  factor  of  about  five  compared 
with  our  earlier  work  by  CC  alone.  This  savings  in  computing  time 
will  increase  progressively  with  number  of  channels  included  in  the  coupled 
equations.  Viewed  from  the  point  of  RM,  it  is  impressive  to  see  that  the 
scattering  functions  can  be  represented  by  as  few  as  two  or  three  basis  functions 
per  channel.  We  believe  this  hybrid  method  will  make  a  very  useful  tool  in 
dealing  with  electron  scattering  from  such  moderate-size  molecules  as  N^  or  C^. 
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Table  I. 

The  number  of  basis  functions  per  channel  used  in  Eq. (6) 
to  calculate  the  A=0  (odd)  partial  cross  section  at  E=15  eV. 

Channel 

£ 

Number  of  basis  functions 

1 

2 

X]Eg 

3 

3 

5 

2 

0 

4 

b  K 

2 

3 

4 

1 

0.000  159 


0.002  684 


Table  IV.  Total  excitation  cross  sections  of  the  b^E,+, 
in  units  of  10'1'  cm2. 
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PART  III 

EFFECT  OF  TARGET  POLARIZATION  IN  SLOW- ELECTRO')  COLLISIONS: 
EXCITATION  CROSS  SECTIONS  OF  H2  (xV  -  bV) 

BY  THE  PO LARI  ZED-ORBITAL  METHOD 
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I.  INTRODUCTION 

In  the  electron- impact  processes  of  atoms  and  molecules,  if  the  energy 
of  the  incident  electron  is  low  enough,  then  the  bound  electrons  of  the 
target  atom  (or  molecule)  adjust  themselves  adiabatically  in  response  to  the 
incoming  electron.  The  distorted  atoms  in  turn  interact  back  with  induced 
multipole  moments  in  form  of  long-range  potentials.*'3 4  This  is  the  character¬ 
istic  feature  of  slow-electron  scattering  in  contrast  to  collision  by  fast 
electrons.  Thus,  simple  theoretical  prescriptions  such  as  the  Bom  approxi¬ 
mation  are  no  longer  valid.  In  principle  the  distortion  of  the  target  may 
be  taken  into  account  by  including  a  sufficient  number  of  terms  in  the  close¬ 
coupling  expansion.  Such  an  expansion  may  be  further  improved  by  including 
pseudostates  in  the  expansion.  This  method  has  shown  to  be  satisfactory 
in  calculating  excitation  cross  sections  of  atomic  H.5 * * 8  Alternatively,  one 
may  direct  attention  to  the  long-range  potentials,  and  construct  empirically 
these  potentials  from  the  known  static  properties  of  the  target  atoms.  The 

latter,  called  pseudopotential  method,  has  been  applied  to  elastic  scattering 
6  7 

problems  on  and  N^. 

3  4 

Another  approach  is  the  method  of  polarized  orbitals,  *  in  which 

adiabatic  first-order  corrections  to  the  target  molecular  orbitals  are 

explicitly  included  as  part  of  the  target  wave  function.  The  polarized  orbitals 

(correction  terms)  are  calculated  with  only  the  dipole  component  of  multi- 

pole  moments  included  as  the  perturbation.  With  the  molecular  orbitals  (MO) 

thus  expressed  as  linear  superposition  of  the  unperturbed  and  polarized 

8 

orbitals,  the  standard  close -coupling  (CC)  calculation  can  be  carried  out 
to  obtain  cross  sections. 

In  this  report  we  calculate  the  excitation  cross  sections  of 
(X*y*  -*  B*j*)  by  a  two-state  close -coupling  with  the  polarized  orbitals 


included.  A  comparison  is  made  with  the  previous  two-state  close-coupling 

8 

calculation  without  polarized  orbitals  to  assess  the  influence  of  the  long- 
range  induced  potentials. 

Sec.  II  contains  a  brief  description  of  the  polarized -orbital  method  and 
numerical  procedure  for  obtaining  the  polarized  orbitals.  In  Sec.  Ill  we 
describe  how  these  polarized  orbitals  are  utilized  to  augment  the  close¬ 
coupling  equations  as  well  as  a  brief  recount  of  the  close-coupling  method. 
The  results  and  their  comparison  with  the  previous  calculation  with  no 
polarization  are  presented  in  Sec.  IV. 
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II.  POLARIZED  ORBITAL 


The  purpose  of  polarized  orbitals  is  to  incorporate  the  characteristic 
feature  of  slow-electron  scattering  into  a  mathematical  formulation,  namely, 
the  distortion  of  target  caused  by  slow  electron  located  at  a  substantially 
far  distance  compared  with  molecular  dimension.  In  this  report  we  consider 
only  the  distortion  of  the  target  molecule  by  the  Coulomb  interaction  of  the 
incident  electron.  Moreover,  since  the  dipole  interaction  dominates  over 
other  higher  multipole  terms,  we  will  keep  the  dipole  term  only.  It  will 
be  seen  below  that  this  approximation  will  also  facilitate  the  solutions 
for  differential  equations  for  the  polarized  orbitals. 

Since  the  theoretical  derivations  leading  to  the  differential  equations 
for  the  polarized  orbitals  are  identical  in  all  essential  detail  to  those 
given  by  Temkin  and  co-workers, ^  a  brief  outline  is  given  below  as  they  relate 

Q 

to  the  present  problem.  Next,  a  numerical  procedure  of  solving  inhomogeneous 
differential  equation  is  described. 


A.  Formulation 

The  Hamiltonian  for  the  system  consisting  of  the  target  molecule  and 
the  scattering  electron  located  at  r^  is 


H  .  =  H  ,  +  H' , 
ad  mol  ’ 


H’  -  -zcIVa''1 


(1) 

(2) 


where  r.,  rD  are  the  position  ve  tors  of  the  two  nuclei,  and  Z  (=1)  is  the 

A  D 

nuclear  charge.  We  construct  a  variational  trial  wave  function  for  the 
ground  state  as 
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*X  =  A  UlOgC^)  +  log'(r1|r3))a  (log(r2)  +  log'  (r2|r3)B] ,  (3) 

where  l°g' (rj I r3)  is  the  first-order  correction  to  the  unperturbed  MO  lo  (r^) 

due  to  the  scattering  electron  located  at  r^.  If  the  energy  of  the  system  is 

minimized  with  respect  to  variations  in  lo  ' ,  the  first-order  term  yields  an 

8 

equation  for  the  polarized  orbital  lo  ' , 

(h  +  V  -  e°)  log '  =  -H'lo  ,  (4) 

with 

h  =  - (1/2) V2  -Z(| r-rA|_1  +  (5) 

v  =  /  la*(r1)|r1-r3r1log(r1)dr1,  (6) 

and  e°  is  the  orbital  energy  of  the  unperturbed  lo  .  We  now  truncate  the 

S 

multipole  expansion  in  Eq.  (2),  retaining  the  dipole  term  only, 


H'  *  (4tt/3)  Z  r.  rl2  Z  Y  *(r.)Y,  fr,) 
^  l  3  q  lq^  V  lqK  3J 


where  is  the  spherical  harmonic.  This  allows  us  to  write  the  polarized 
orbital  log(r.jJr3)  as  a  product  of  functions  of  r^  and  of  r^,  viz., 

Substitution  of  Eqs.  (7)  and  (8)  into  Eq.  (4)  yields  a  set  of  coupled 
differential  equations  in  0^  for  each  m. 

(~2  -  ^  ■  h'  Vfr>V.<r>  *  V)'  (« 


UM,M  ■  -2Z  /df  Ohtlr-^r1  *  |r-r„|  'hY-.fi) 


+  2  /d?’dr  loe(r')Y  (rljr-^'l'ho  (?’)Y  ,  (r), 


(10) 
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Sn(r)  =  -(8^3)r  /d?,lag(?,)Yi,-ra^')Y£m^')»  (ID 

and 

k2=2|e°|.  (12) 

The  differential  equations  in  Eq.  (9)  resemble  close-coupling  equations  except 
for  the  inhomogeneity  G^Cr) ,  which  represent  the  dipole  component  of  the 
perturbation  caused  by  the  incident  electron.  The  numerical  procedure  used 
for  the  solutions  of  these  equations  is  described  in  the  next  subsection. 


B.  Numerical  Solution 
Let  us  define 

\v  -hv  *  [^-k2]«u..  <13> 

and  rewrite  Eq.  (9)  in  a  matrix  form 
A2 

— $  =  W<J>  +  G,  (14) 

dr  "  ~ 

where  $  and  G  are  viewed  as  column  vectors  and  W  a  matrix.  The  physically 
acceptable  solutions  $  must  satisfy  the  boundary  conditions 


* 

o 

II 

/- - N 

o 

V _ ^ 

■©- 

(15) 

-Vj* 

<^(r)  ^  e  as  r-**>. 

(16) 

The  numerical  method  employed  in  solving  Eq.  (14)  is  similar  to  the  one  used 
by  Vo  Ky  Lan,^  and  consists  of  replacing  the  differential  equations  by  the 
difference  equations,  i.e., 


ai  ti.k  ■  *  5k  ■ 


wtere  ^  denotes  (Hr^)  etc.,  and  a1)  is  the  Lagrange  n-point  differentiation 
coefficient  at  the  point  r^.  Next,  we  assume 
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sj  =  o,  Cq  =  o.  (22) 

Starting  from  k  =  0,  the  recursion  relations  of  Eqs.  (20)  and  (21)  are  used  to 

generate  arrays  and  ^  to  a  suitably  large  r^_^ast .  Then  we  set  4>^ast  =  0, 

and  the  inward  recursion  formula  of  Eq.  (18)  is  used  to  generate  all 

In  our  work  we  adopted  five-point  Lagrange  formula  (n=5)  in  Eq.  (17)  as 

9 

opposed  to  the  three-point  formula  used  by  Vo  Ky  Lan.  But  there  is  no 
significant  difference  between  the  two  schemes.  The  outward  recursion  was 
carried  out  to  23  -  33  a  ,  Satisfactory  results  are  obtained  with  r,  in  this 

O  n 


range. 


III.  SCATTERING  EQUATIONS 
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The  scattering  equations  for  a  close  coupling  is  well  known  in  the 
case  where  unperturbed  MO  are  used  only,  i.e., 


2 

[d  _  +  k2]  F  ,  (r)  =  E  „(U  „  +  W  „)  F  „  ,  (r) , 

.  Z  Z  pJ  PM  ’  ii"  mi"  mi"''  n"n,v  ’ 


p"  pp' 


PP"'  p"p' 


where  U  „  and  W  „  are  direct  and  exchange  potentials  respectively.  If  the 
polarized  orbitals  are  included  in  the  target  state  in  the  manner  of  Eq.  (3), 
then  additional  polarization  potentials  enter  Eq.  (23)  in  an  analogous 
way  as  U  „  and  W  ... 

7  pp"  pp" 

However,  before  plunging  into  calculation,  it  is  worthwhile  to  review 
what  we  are  attempting  to  accomplish  with  the  polarized  orbitals  in  the 
context  of  physical  situation.  Polarized  orbitals  are  to  describe  the 
distortion  of  target  molecule  when  the  incident  electron  is  at  a  distance  sub¬ 
stantially  greater  than  the  molecular  dimension.  When  the  incident  electron 
comes  within  the  molecule,  all  electrons  including  the  incident  electron  must 
be  treated  on  an  equal -footing.  This  physical  situation  cannot  be  described 
by  polarized  orbitals.  Therefore,  it  has  been  tacitly  assumed  that  the 
the  incident  electron  always  stays  outside  the  molecular  electrons.  We  now 
make  this  assumption  explicit  as  suggested  in  Temkin's  original  paper. 

For  example,  in  the  expression  of  the  direct  potential  due  to  the  polarized 
orbital, 


Or)  *  2  El*(,’/2l4l>  /d;Y*,A-A*YL«WYll'A-X'WYln,W 

r 

x  [r"L1  /  d?log(?’)r’Llo-  (r' )Y^(r') 

+  rL  /  drlOgCrnr'^'hog'Cr'JY^)],  (24) 


the  second  integral  is  emitted  because  it  is  the  contribution  when  the  incident 
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k 


electron  lies  inside  the  molecular  electrons.  On  a  similar  ground  the  exchange 
interactions  are  ignored  from  consideration.  Furthermore,  we  have  omitted 
the  coupling  potentials  between  the  ground  and  excited  states  due  to  polarized 
orbitals.  This  is  consistent  with  the  physical  description  of  the  incident 
electron  being  at  a  far  distance  where  it  "sees"  the  molecule  in  one  or  the  other 
state. 

The  procedures  of  solving  the  CC  integro- differential  equations  and  of 
extracting  cross  sections  from  their  asymptotic  form  are  described  in  detail 

g 

elsewhere,  and  will  not  be  repeated  here. 


IV.  RESULTS  AND  DISCUSSION 


SI 


The  methods  described  in  Secs.  II  and  III  are  used  to  calculate  the 
excitation  B2E^)  cross  sections  of  molecule.  The  configurations 

of  these  states  are: 


(la  )2 

g 

(!ag)  (lauJ . 


The  unperturbed  MO's  are  obtained  by  two  separate  SCF  calculations  and  identical 

8  + 

to  those  used  in  our  previous  work.  Polarized  orbitals  of  o  ,  tt  .  it 

u’  u’  u 

symmetries  are  calculated  from  the  parent  la  of  the  ground  state,  and  those 

g 

of  a  ,  tr  ,  it  symmetries  are  calculated  from  la  orbital.  Polarized  orbitals 

o  o  o  ^ 

from  the  parent  la  of  the  B^E*  are  also  calculated.  But,  since  they  are  much 

o  u 

smaller  in  magnitude  than  those  from  the  la^  orbital,  we  have  excluded  them 
from  further  consideration.  This  is  understandable  because  the  la 


(e  =  -24.94  eV)  orbital  of  the  B1!*  state  is  much  more  tightly  bound  than 
lo  (e  =  -5.84  eV)  and  expected  to  be  less  susceptible  to  polarization.  We 
obtained  4.78  a.u.  for  the  isotropic  component  of  the  dipole  polarizability 
of  Wj,  which  is  within  8%  of  the  best  calculated  value  5.18  a.u.  available.^1 
This  good  agreement  gives  an  independent  check  for  the  accuracy  of  our  polarized 
orbitals. 

The  excitation  cross  sections  are  obtained  at  two  representative  energies 
of  the  incident  electron  and  shown  in  Table  7.  For  comparison  we  included  in 
Table  I  cross  sections  obtained  without  polarized  orbitals.  The  inclusion  of 
polarized  orbitals  makes  a  significant  (>35%)  difference  in  the  cross  section 
at  E  =  25  eV,  but  much  less  change  (9%)  at  E  =  50  eV.  This  is  consistent 
with  the  physical  reasoning  that  polarization  is  a  phenomenon  of  slow- 
electron  scattering. 


► 


L 
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Unfortunately,  the  only  available  set  of  experimental  cross  sections^ 
are  in  rather  substantial  disagreement  from  the  close -coupling  calculation 
with  or  without  polarized  orbitals  included.  One  can  only  hope  that  independent 
measurements  and  calculations  will  clarify  the  situation. 

From  this  study  it  is  apparent  that  the  target  polarization  is  an 
important  factor  to  be  considered  in  the  collision  processes  of  molecules  by 
slow  electrons.  However,  theoretical  means  to  achieve  this  goal  needs  further 
improvement . 
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Table  I.  Direct  excitation  cross  sections  of  the  state  of  H 
in  units  of  10-17  cm2.  u 


Incident 

Close  Coupling  with 

Standard 

Experiment' 

Energy  (eV) 

Polarized  Orbitals 

Close  Coupling 

25 

5.8 

4.3 

2.0  ±  0. 6b 

50 

5.1 

4.7 

2.8  ±  0.8 

aRef .  11 

^Linear  interpolation  of  measured  values  1.9  ±  0.6  at  20  eV  and  2.0  ±  0.6 
at  30  eV. 
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PART  IV 

TRANSITION  PROBABILITY  OF  THE  1SQ  -  ^ 
EMISSION  OF  ATOMIC  OXYGEN 
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I.  INTRODUCTION 


The  dipole- forbidden  transition  -»■  ^  (5577  8)  in  01  is  a  prominent 
feature  of  auroral  spectra,  which  has  received  a  good  deal  of  attention 
historically.  The  ->  ^  transition  probability  also  plays  an  important 
role  as  a  parameter  in  quantitative  modeling  of  the  earth's  upper  atmosphere, 
especially  in  establishing  population  densities  of  various  atomic  and  molecular 

species  and  in  determining  mechanisms  of  production  and  quenching  of  them.1 

1  3 

The  upper  state  S„  also  makes  a  magnetic  dipole  transition  to  the 

state  (A  =  2972  8)  besides  the  5577  8  emission.  The  transition  probabilities  of 

these  lines  as  well  as  the  lifetime  of  the  ^  state  have  been  determined 

2_  s 

experimentally  by  numerous  groups.  The  experimental  value  of  the  Einstein  A 

coefficient  for  the  5577  8  line  is  1.06  ±  0.32  sec  \  and  the  decay  constant 

of  the  1S()  state  [i.e.,  A(5577  8)  +  A(2972  8)1  ranges  from  1.11  to  1.43  sec'1.3 

9-13 

There  have  also  been  several  theoretical  calculations  with  the  resulting 
value  of  A  between  1.183  and  2.2  sec  1.  The  earlier  calculations  in  which 
hydrogenic  wave  functions  were  employed  are  thought  to  have  only  limited 
accuracy. More  recently,  Nicholaides  et  al.13  obtained  A(5577  8)  as 
1.183  sec  1  based  on  wave  functions  computed  by  the  method  of  "nonclosed- 
shell  many-electron  theory'  (NCMET) ." 

'The  effects  of  electron  correlation,  which  were  important  distinguishing 

features  of  the  work  by  Nicholaides  et  al_ . ,  may  be  taken  into  account  in  a 

14 

systematic  way  by  the  multiconfiguration  Hartree-Fock  (MCHF)  method.  In 
view  of  the  importance  of  the  -+  ^  transition  in  atmospheric  research, 
it  is  worthwhile  to  reexamine  the  theoretical  transition  probability  from 
the  MCHF  approach. 


II.  METHOD  OF  CALCULATION  AND  RESULTS 
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The  wave  functions  employed  in  this  work  have  been  computed  by  the  MG  IF 
14 

method,  and  are  expressed  in  terms  of  linear  combinations  of  single  con¬ 
figuration  functions  </k  as 


n^) 


E i  ai4/i(r1,r2,....,i 


V » 


(1) 


and 

v(1D2)  =  Sibi^i(r1,r2,....,rNi1D2),  (2) 

where  a^,  b^  are  configuration -mixing  coefficients.  The  MCHF  method  is 
described  in  detail  elsewhere.14  Two  sets  of  MCHF  functions  of  the 
state  are  displayed  in  Table  I  showing  the  configurations  and  the  corresponding 
coefficients  as  well  as  the  computed  total  energy  for  each  of  the  two  sets. 

In  the  cases  where  more  than  one  LS  coupling  term  may  be  derived  from  an  (nd)^ 
configuration,  a  specific  designation  is  written  inside  a  square  bracket. 
Similarly,  we  show  two  sets  of  the  1D2  functions  in  Table  II.  In  the  present 
MG1F  procedure,  the  orbitals  and  the  configuration  coefficients  are  optimized 
simultaneously  so  that  the  orbitals  are  common  for  all  configurations  within 
a  state.  Although  no  relationship  need  be  imposed  between  sets  of  orbitals 
belonging  to  one  state  and  of  those  belonging  to  the  other  state,  the  Is,  2s, 
and  2p  orbitals  are  held  common  to  both  the  ^  and  1D.,  states  so  that  we  may 
take  advantage  of  the  orthonormal  property  of  the  orbitals  when  evaluating 
the  quadrupole  matrix  elements  between  the  1Sq  and  4D2  states.  Even  with  this 
restriction,  we  believe  the  wave  functions  are  quite  accurate  as  a  sufficient 
number  of  configurations  are  included  to  account  for  effects  such  as  electron 
correlations.  For  example,  the  computed  1SQ  -  JD2  energy  splitting  is 
0.086465  a.u.  (5270  8)  from  set  A  of  Tables  I  and  II,  which  is  within  5.5".  of 
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the  observed  value  of  5577  8.  In  studying  the  forbidden  transitions,  the 
fundamental  quantity  is  the  quadrupole  line  strength  S  defined^  as 


S(1S(),  !D2)  =  (2/3je2  Em|(W5)1/2//(?1,?2,..JN|]S0); 


r's  rs  Y2m(V  '«'(r1,r,>..,rM|1D,)dr1dr,...drk;|2, 


12' 


(3) 


where  is  the  sperhical  harmonic.  The  summation  over  m  is  merely  formal 
since  there  is  only  one  nonvanishing  element,  and  the  summation  over  s  covers 
all  eight  atomic  oxygen  electrons.  In  the  LS- coupling  scheme,  the  quardupole 
transition  probability  is  defined  through  the  line  strength  as^’^ 

A(1Sq  +  1D2)  =  327t6S(1S0,  ^J/ShX5.  (4) 

17 

A  related  quantity  is  the  quadrupole  transition  oscillator  strength  defined 
as 

f  (1S0,1D2)  =  4Tr4mcS(1S0,1D2)/5he2A3,  (5) 


where  m  and  c  are  the  mass  of  electron  and  the  speed  of  light  respectively. 
With  the  wave  functions  in  the  MCHF  form  as  in  Eqs.  (1)  and  (2),  the 
line  strength  is  now  written  as 


SC1S0,1D2)  =  (2/3)e2|EijIij|2,  (6) 

with 

1  ij  =  aibj(47T/5)1/2/^*(r1,r2,...,rN|1S0)x 

Es  rsY2m(rs)Vrl,,‘,rNl  D2)dr1dr2...drN,  (7) 

where  m  is  to  take  the  particular  value  (from  -2  to  2)  that  gives  a  nonvanishing 
integral.  In  Table  III  we  display  the  contributions  h j  computed  by  Eq.  (7) 


mil  if liiVim  aaifruMfciili 
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along  with  the  wave  functions  of  set  A  in  Tables  I  and  II.  We  note  from 

Table  III  that  the  primary  configuration  ls22s22p4  alone  gives  2.51  a2  as 

2 

compared  with  2.42  aQ  when  all  contributions  are  summed  (a^  being  the  Bohr 

radius).  The  corresponding  values  with  the  wave  functions  of  set  B  are 
2  2 

2.55  aQ  and  2.48  aQ  respectively.  The  good  agreement  (31)  between  the 

results  from  the  larger  (A)  and  smaller  set  (B)  sets  indicates  that  the  wave 

functions  of  either  set  are  reasonably  accurate.  By  using  our  best  value 

(2.42  a2)  and  with  A  =  5577  ft,  we  obtain  1.22  sec  1  for  the  transition 

-9 

probability  and  5.69  x  10  for  the  oscillator  strength. 

From  the  off-diagonal  hypervirial  relation,^ 

£s  v2s  =  (2/AE)  Es  rs-?s  .  (8) 

we  obtain  an  equivalent  expression  of  Eq.  (7),  i.e., 

lij  =  ai.bj  (4t/5)1/2  (2/AE)  /  (r^, . .  ,rN|  \)* 

ls  Y2jiSTs^  rs'^s  '*Jj<'rl,r2’*",rN^  °2^  drldr2 •  •  •  drN » 

where  AF,  is  the  energy  difference  in  atomic  units.  When  Eq.  (9)  is  used  to 

compute  the  integrals,  the  resulting  strength  is  in  the  "velocity  form"  as 

opposed  to  the  "length  form"  when  Eq.  (7)  is  used.  We  re-compute  the 

oscillator  strength  with  the  velocity  form  and  include  in  Table  III  the 

contributions  from  the  various  configurations.  The  velocity  form  gives  an 

-9 

oscillator  strength  of  4.88  x  10  in  reasonable  agreement  with  the  result 

of  the  length  form.  Furthermore,  in  the  velocity  form  there  is  no  contribution 

2  2  4 

to  the  line  strength  from  the  dominant  Is  2s  2p  configuration  (Table  III)  so 
that  this  agreement  is  achieved  through  the  inclusion  of  correlation.  Thus 
the  velocity  form  provides  a  critical  test  of  the  accuracy  of  the  wave 
function.  Both  these  observations  render  additional  support  to  the  accuracy 
of  the  wave  functions  used  in  this  work. 


III.  DISCUSSION 


The  transition  probability  (1.22  sec  *)  of  this  work  agrees  very 
well  with  1.28  sec  ^  computed  by  Garstang^’^  and  1.183  sec  ^  by 
Nicholaides  et  ad.^  Garstang's  value  of  1.28  sec  ^  is  based  on  the  single¬ 
configuration  Hartree-Fock  functions.  This  is  in  general  accord  with  our 

2  2  4 

finding  that  in  the  length  form  the  primary  configuration  of  Is  2s  2p 

dominates  over  all  others  as  shown  in  Table  III.  The  close  agreement  with 

the  result  of  Nicholaides  et  al.  is  particularly  encouraging,  since  their 

method  of  computing  the  wave  functions  (NCMET)  is  quite  different  from  the 

12 

MCHF  method  adopted  in  this  work.  In  Garstang's  semiempirical  study,  the 
effect  due  to  the  departure  from  the  LS-coupling  scheme  is  found  to  be  very 
small;  1.25  sec  1  and  1.28  sec  1  respectively  with  and  without  such  an  effect. 
Therefore,  it  is  unlikely  that  the  theoretical  value  presented  here  will 
change  significantly  even  with  the  inclusion  of  other  effects  such  as  spin- 
orbit  coupling. 

The  experimental  determination  of  the  transition  probability  A(5577X) 

is  impeded  by  the  difficulty  of  ascertaining  the  number  density  of  the 

state  and  by  the  competing  2972$  (  -*•  transition.  However,  McConkey 

and  kemahan^  measure  A(5577^)  to  be  1.0  sec  *  from  a  discharge  experiment. 

By  using  a  similar  technique,  Kemahan  and  Pang'  obtained  a  value  of 

1.06  f  0.32  sec  *.  The  greater  certainty  in  the  latter  experiment  is 

largely  due  to  a  more  precise  knowledge  of  the  number  density.  Other 

experimental  data  relevant  to  this  work  are  measurements  of  the  radiative 

decay  constant  of  the  state,  i.e.,  the  sun  of  A(S577$)  and  A(2972^) . 

8  -I 

Gomey  and  Williams  obtained  1.31  ±  0.05  sec  for  the  decay  constant  T  from 
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a  laboratory  experiment.  From  the  analysis  of  aurora  observation,  Omholt4, 
and  Evans  and  Jones‘S  deduced  F  as  1.43  and  1.49  sec  *  respectively.  In 

addition  the  A(5577$)/A(2972$)  ratio  has  been  reported  as  18.6  by  McConkey 

->  6  2 
et  al . '  ,  22  by  Le  Blanc  et  al.  ,  and  23.7  by  Kemahan  and  Pang.  From  the 

experimental  data  of  the  decay  constant  and  of  the  A(5577S)/A(2972$)  ratio, 

we  estimate  the  value  of  ACSS???''  to  be  in  the  range  of  1.24  sec  1  to  1.37  sec  * 

-1  2 

The  direct  measurement  1.06  ±  0.32  sec  by  Kemahan  and  Pang  is  some¬ 
what  smaller  than  our  theoretical  value  1.22  sec  though  the  latter  is  well 
within  the  stated  uncertainty  of  the  experimental  data.  On  the  other  hand 
the  "deduced  values"  (1.24  -  1.37  sec  *)  are  somewhat  greater  than  the 
theoretical  value.  Considering  the  experimental  difficulties  and  uncertainties, 

it  is  encouraging  to  see  this  essential  agreement  between  theory  and  experiment. 

11-13 

The  good  agreement  of  our  calculation  with  the  previous  theoretical  works  , 
where  realistic  wave  functions  were  used,  also  gives  a  greater  confidence  in 
the  A(55778)  value  presented  here. 


The  research  presented  in  Part  IV  was  done  in  collaboration  with 

Professor  Charlotte  F.  Fischer  of  the  Vanderbilt  University.  The  part 

of  the  work  contributed  by  Professor  Fischer  was  not  supported  by  the 

Air  Force  Geophysics  Laboratory.  A  paper  entitled  "Transition  Probability 
1  1 

of  the  SQ  -i-  Emission  of  Atonic  Oxygen"  has  been  published  jointly 
with  Professor  Fischer  in  the  Journal  of  Chemical  Physics,  Volume  76, 
pp.  498-501(1982). 
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Table  I.  Mult iconfigurat ion  Wave  Function  of  Oxygen  *S  State. 


Configuration 

Coefficients 

Set  A a)  b) c) 

Set 

R  c) 

1. 

(2s)2(2p)4 

0.970  206 

0.978 

474 

2. 

(2p)  6 

0.189  263 

0.201 

885 

3. 

(2p)4[lD](3Pl)2 

0.005  880 

0.010 

334 

4. 

(2s)2(2p)3[2P](3p2) 

0.018  115 

-0.011 

121 

5. 

(2p)5(3p3) 

-0.028  477 

0.030 

984 

6. 

(2s)(2p)4[1D](3d1) 

0.030  612 

0.025 

077 

7. 

(2s)(2p)4[1D](3s1) 

0.000  625 

0.003 

363 

8. 

(2s)2(2p)2[1S](3p4)2[1S] 

-0.072  190 

9. 

(2s)2C2p)2[3PJ(3par  [3P] 

0.053  232 

10. 

(2s)2(2p)2[1Dj(3p4)2[1D] 

0.059  830 

11. 

(2s)2r2p)2[1S](3s2)2 

0.022  216 

12. 

f2p)4[1S](3s2)2 

-0.019  487 

13. 

f2s)2(2p)2[1S](3d2)2[1S] 

0.067  818 

14. 

(2s)2(2p)2[5P](3d2)2[3P) 

-0.044  863 

13. 

(2s)2(2p)2[1Dj(3d2)2[1D] 

-0.040  721 

Hnergy 

-74.759  539  a.u. 

-74.670 

000  a 

a)  Subscripts  indicate  different  orbitals  of  the  same  type. 

b)  Is,  2s,  and  2p  orbitals  from  fully  variational  MOIF  calculation  for  first  three 
configurations;  all  other  orbitals  variational. 

c)  Is,  2s,  2p  and  3p^  orbitals  from  fully  variational  MCltF  calculation  for  first 
three  configurations;  all  other  orbitals  variational. 


Table  11.  Multiconfiguration  Wave  Function  of  Oxygen  1)  State 


Configurat ion 


a) 


Coefficients 


Set  A 

Set 

B  c) 

1. 

(2si2(2p)4 

0.984  825 

0.993 

420 

2. 

Up)  Vb]  <3Pj)2 

0.01b  950 

0.016 

257 

3. 

(2s)\2p)3[2P](3p2) 

-0.020  586 

-0.028 

079 

4. 

(2s)2(2pl5l2D|(3p2) 

-0.000  783 

-0.006 

358 

5. 

(2p)5(3p3) 

-0.029  366 

-0.029 

958 

6 . 

(2s)(2p)4l1S](3d1) 

0.000  962 

-0.000 

508 

7. 

(2s)(2p)4f3P](3di) 

-0.103  608 

-0.105 

378 

8. 

(2s)i2p)4liD](3d1) 

0.002  907 

-0.000 

497 

9. 

(2sj(2p]4fVl(3s1j 

-0.005  460 

-0.004 

695 

10. 

(2s)2(2p)2[1S](3p)2[1D] 

0.02  6  074 

11. 

(2s) 2  (2p) ^  [3P]  ( 3p 4 )  2  f^P] 

0.054  850 

12. 

(2s) 2 (2p) 2 [ lDI  (3p4) 2  f 3P] 

C. 035  520 

13. 

(2s)2(2p)2!ID](3p4)2[1l)] 

0.070  313 

14. 

(2s) 2  (2p)  2  ( 1 D ]  (3s  2) 2 

-0.010  932 

15. 

(2p)4 [ JD] (3s2)2 

-0.023  824 

16. 

(2s)2(2p)2[1S](3d2)2f1D] 

-0.015  317 

17. 

(2s)2(2p)2[3Pl (3d2)2[3P] 

-0.045  117 

18. 

(2s.)2(2p)2[1D](3d2)2i1S] 

-0.034  052 

19. 

(2s)2(2p)2[1D](3d2)2f1I)] 

0.043  443 

20. 

(2s)2(2p)3[2P](4f) 

0.024  741 

21. 

(2s) 2 (2p) 3 ( 2  D 1 (4f) 

0.035  528 

Energy 

-74.846  004  a.u. 

-74.770 

167  a.u. 

a)  Subscripts  indicate  different  orbitals  of  the  same  type. 

b)  Is,  2s,  2p  orbitals  fixed  as  in  Set  A  of  Table  1;  all  other  orbitals  variational. 

c)  A&,  2s,  2p  and  3p^  orbitals  as  in  Set  3  of  Table  1;  all  other  orbitals  variational 


Table  III. 

Contributions  to  quadrupole  matrix  element  I 

65 

Configurations 

1  . 

1.1 

lS 

lD2 

by  Eq. 

(7) 

by  Eq. 

(9) 

1 

1 

2.511 

346 

0.0 

1 

3 

-0.062 

992 

-0.528 

617 

1 

6 

-0.002 

412 

-0.054 

722 

1 

20 

-0.092 

596 

-2.381 

701 

2 

S 

0.020 

423 

0.163 

222 

4 

1 

0.026 

269 

-0.162 

122 

4 

3 

-0.002 

158 

0.004 

720 

4 

20 

-0.001 

001 

0.007 

322 

5 

5 

-0.002 

600 

0.000 

022 

6 

1 

-0.029 

648 

0.741 

202 

8 

3 

0.002 

121 

-0.021 

689 

8 

10 

0.009 

408 

-0.001 

819 

8 

13 

0.006 

737 

0.0 

9 

3 

-0.001 

173 

0.011 

995 

9 

11 

0.011 

132 

-0.001 

410 

10 

10 

-0.001 

833 

0.0 

10 

12 

-0.004 

750 

0.000 

918 

10 

13 

0.018 

978 

-0.002 

404 

12 

15 

0.001 

220 

0.0 

13 

16 

0.006 

444 

-0.001 

246 

13 

18 

0.006 

069 

0.0 

14 

17 

0.010 

087 

-0.001 

435 

IS 

18 

-0.003 

847 

0.000 

744 

Others 

-0.000 

282 

-0.018 

289 

TOTAL 


2.424  942 


-2.245  309 
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